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in every energy range and the evolution of /Var with the energy shows a similar increase both in 
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1. Introduction 

As many blazars, PKS 2155-304 (z = 0.117) features a questioning and interesting behavior: 
its variability which can be seen over a great energy range, from radio to TeV. Characterizing 
the variability of blazars is important to study the underlying acceleration processes which are, 
at the moment, poorly known. Multiwavelength (MWL) data, from experiments and instruments 
as H.E.S.S., Fermi-LPLL, RXTE, Swift-XRT, XMM-Newton and SMARTS are used to study the 
variability of this object. 

Lognormality and intrinsic variability of the MWL light curves are investigated. Lognormal 
flux variability was first found in galactic accreting sources like X-ray binaries, linking lognormal 
processes to accretion disks [1]. In a lognormal process, the fluctuations of the flux are on average 
proportional, or at least correlated, to the flux itself, ruling out additive processes in favor of multi¬ 
plicative ones. Lor blazars, a lognormal flux behaviour could as well be indicative of the accretion 
disk variability’s imprint onto the jet [2]. BL Lacertae was the first blazar to show a lognormal be¬ 
havior in the X-ray [3]. Lor PKS 2155-304, the first detection of a lognormal behavior was during 
the flaring state of 2006 in H.L.S.S. data [4]. 

Studying the evolution of the variability with the energy is one element to further distinguish 
between different acceleration models. 

2. Datasets 

2.1 Very Bight Energy (VHE): H.E.S.S. 

The H.L.S.S. (High Energy Stereoscopic System) experiment is an array of 5 ground-based 
Cherenkov Telescopes which detect the Cherenkov light emitted by the shower created by VHE 
y-rays with energies from ^100 GeV to > 10 TeV. The four first telescopes are operated since 
2003 (HESS-I phase) and a fifth telescope had been added to the array in 2012 (HESS-II phase). 

The present analysis is based on the data taken between MJD 53200 and 56246 using only 
the HESS-I phase of the array (without the fifth telescope) collected with 3 or 4 telescopes. The 
flaring state of July/August 2006 from MJD 53943 to 53949 has been removed, in order to have a 
light curve of the quiescent state of PKS 2155-304 [4]. In total, about 260h of data passed standard 
quality cuts as defined in [5], with a mean zenith angle of 21° and an energy threshold of ISOGeV. 
The data are analyzed with the Model analysis chain (Model-i-i-) [6] between 0.2 and 125 TeV, 
yielding to a total detection significance >300(7. 

The light curve of nightly fluxes used here is calculated assuming a power-law spectrum 
dN/dE oc E~^, r being the power law index. However, it has been reported in [4] that PKS 2155- 
304 shows some spectral variability. Therefore the fluxes are derived with a different fitted spectral 
index each year to take into account this variability. This leads to 9 different values for L (one per 
year) with a mean of 3.45 zb 0.05 and a standard deviation of 3%. 

The resulting quiescent nightly light curve, presented in the first panel of figure 1, has a mean 
flux of 6.9 X 10“^^ph.cm“^.s“^ This analysis procedure has been applied to the same data set 
using the following energy cuts: 0.2-0.63 and 0.63-125 TeV, the cut at 0.63 TeV being chosen to 
have the same statistics in each energy band. 
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Figure 1: Light curves from top to bottom : H.E.S.S. in cm“^.s“^ from 0.2 to 125 TeV, Fermi-LP4L in 
cm“^.s“^ from 0.1 to 300 GeV, X-ray in erg.cm“^.s“^ from 2 to 10 keV (Swift-XRT in black, RXTE in blue 
and XMM-Newton in green) and SMARTS (R band) in Jy. 


2.2 High Energy (HE): FermEhAT 

The Fermi Large Area Telescope (LAT) [7] on board the Fermi satellite is a pair-conversion 
telescope designed to detect HE y-rays in the energy range from ^ 20 MeV to more than 500 GeV. 

The data for PKS 2155-304 used here goes from MJD 54687 to 56687 and are taken be¬ 
tween 100 MeV and 300 GeV. The detector is described by the P7REP_SOURCE_V15 instrument 
response functions^ and the analysis has been performed by the tool enrico [8] using Fermi Sci¬ 
ence Tools v9r32p5. The prefactor of the diffuse galactic background and the normalisation of the 
isotropic diffuse emission are left free to vary in the likelihood fit. All sources from the 3FGL 
catalog [9] within 15° around PKS 2155-304 are included in the model. The spectra of PKS 2155- 
304 and the nearby FSRQ (Flat Spectrum Radio Quasar) PKS 2149-306 are modeled by a simple 
power-law, and their spectral indices and prefactors are left free. All other components are fixed to 
the values in the 3FGL catalog. 

PKS 2155-304 is detected with a significance of > 150(7 with an overall spectral index of 

^http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/Cicerone_Data_Exploration/ 
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1.83 ±0.01. The photon counts are integrated into 201 bins of 10 days, resulting in a light curve, 
presented in the second panel of figure 1, with a mean flux of 1.3 x 10“^ph.cm“^.s“^ The same 
analysis procedure has been applied to the same data set using the following energy cuts : 0.1-3, 
3-30 and 30-300 GeV. 

2.3 X-ray: RXTE^ Swift-XRT and XMM-Newton 

RXTE (Rossi X-ray Timing Explorer) already analyzed light curves are publicly available^ 
between 2 and 10 keV, and in three other energy bands: 2-4, 4-7 and 7-10 keV. 

Swift X-ray Telescope (XRT) data were analyzed using the package HE AS OF T 6.16 between 
2 and 10 keV, and also in four other energy bands: 0.3-2, 2-4,4-7 and 7-10 keV. The data were recal¬ 
ibrated using the last update of CALEB and reduced using the standard procedures xrtpipeline 
and xrtproducts. Caution has been taken to properly account for pile-up effects for corre¬ 
sponding affected exposures, and spectral fits were performed using Xspec 12.8.2 assuming a 
power-law spectrum. 

XMM-Newton (X-ray Multi Mirror Mission) public EPIC (European Photon Imaging Cam¬ 
era) data were reduced using the SAS software package, version 14.0, and analyzed following [10] 
between 2 and 10 keV, and also in four other energy bands: 0.3-2, 2-4, 4-7 and 7-10 keV, to match 
RXTE and Swift-XRT data. 

Swift-XRT and XMM-Newton data have been corrected from galactic absorption with nu = 
1.48X [11], 

In all common energy bins, data of the three instruments are combined to construct an all-over 
X-ray light curve. In order to have an X-ray light curve of the quiescent state of PKS 2155-304, 
all ToO (Target of Opportunities) data were removed. The final light curve, presented in the third 
panel of figure 1, goes from MJD 53501 to 55474 with a mean flux of 2.7 x 10“^^erg.cm“^.s“^ 

2.4 Optical wavelengths : SMARTS 

SMARTS (Small and Moderate Aperture Research Telescope System) already analyzed data 
are publicly available^. Magnitudes are corrected for the absorption of the galactic foreground 
using [12] before being converted to a spectral flux density using zero magnitude flux values of 
[13]. The light curves are taken in the J, R, V and B bands from MJD 54603 to 56622. Eor the sake 
of simplicity, only the light curve for the R band is presented in the last panel of figure 1. 

3. Results 

3.1 Lognormality 

The lognormal behavior of PKS 2155-304 is investigated by looking at the distribution of 
fluxes and the correlation between the excess rms and the average flux. 

The flux and log flux distributions of each light curve are fitted by a Gaussian using a x'^ fib 
The results of the fits are summarized in table 1. In all cases, a Gaussian fits the log flux distribution 


^http://cass.ucsd.edu/ rxteagn/ 

^ http ://w w w. astro .y ale. edu/smarts/glast/home.php# 
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;i^^/d.o.f. 

Prob 

log O 

;i^^/d.o.f. Prob 

H.E.S.S. (> 200GeV) 

2.71 

10-4 

0.716 

0.75 

Fermi-LAT (0.1-300GeV) 

2.51 

10-3 

1.36 

0.18 

X-ray (2-lOkeV) 

6.67 

10-10 

1.36 

0.20 

SMARTS (B band) 

4.35 

10-3 

0.422 

0.96 

SMARTS (V band) 

2.30 

10-*^ 

0.790 

0.65 

SMARTS (R band) 

5.42 

10-9 

0.823 

0.62 

SMARTS (J band) 

2.33 

10-3 

1.28 

0.22 


Table 1: Values of the and associated probability (probability that the model is good) for the Gaussian fit 

of the fiux and log fiux distributions for each light curve. 


better than the flux distribution, indicating that MLW fluxes of PKS 2155-304 follow a lognormal 
distribution. 

The scatter plots of the excess rms Oxs = — o^, being the variance and the mean 

square error of the light curve [15], versus the mean flux ^ are shown for each light curve in figure 
2. Each point is calculated with at least 20 light curve points. To test the rms-flux relation, scatter 
plots are fitted by a constant and a linear adjustments. The fit results are summarized in table 2. The 
results show more than a 5(7 preference for the linear fit except for the Fermi-hAI data. However 
the correlation between <7xs and ^ is not expected to be linear, especially when the points are 
spread in the plot, so the correlation factor p and the Kendall rank T (which measure the ordering 
of the points [14]) are calculated. In all cases p >0.80 and T > 0.65 meaning that (Jxs and ^ show 
a strong correlation, leading to the conclusion that the fluctuations of the flux are correlated and/or 
proportional to the flux. 



cst 

cst+a X X 

a 

P 

T 

H.E.S.S. (> 200GeV) 

4.15 

0.773 

5.25 

0.84±0.15 

0.72 ±0.27 

Fermi-FKT (0.1-300GeV) 

0.937 

0.107 

2.75 

0.93±0.19 

0.69 ±0.25 

X-ray (2-lOkeV) 

78.2 

28.9 

23 

0.85 ±0.02 

0.67 ±0.24 

SMARTS (B band) 

4222 

892 

225 

0.92 ±0.03 

0.68±0.18 

SMARTS (V band) 

4472 

810 

236 

0.92 ±0.02 

0.75 ±0.19 

SMARTS (R band) 

7542 

946 

326 

0.93 ±0.02 

0.74±0.18 

SMARTS (J band) 

3339 

398 

197 

0.94 ±0.03 

0.80±0.21 


Table 2: Values of the reduced of the constant and linear fits of the scatter plots shown in figure 2 for 
each light curve, with values for the significance (7, the correlation factor p and the Kendall rank T. 


3.2 Variability energy distribution (VED) 

The fractional variability Fyar is used to quantify the variability of the flux taking into account 
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Figure 2: Scatter plot of the excess rms and the average flux for the MWL data. A linear fit is shown in red. 


the errors [15]: 


/Var = 


1 52 - 




± 


A 




2 


It is computed for each data set at the different energies presented in section 2. The VED plot 
is shown on figure 3. This figures shows that the SMARTS variability level is somewhat constant 
while a clear increase is visible in the X-ray band. Then the variability is smaller in the Fermi-LAT 
range and increases again in the H.E.S.S. range. 


4. Conclusion 

The distribution of the flux of the quiescent state of PKS 2155-304 follows a lognormal dis¬ 
tribution and the variations of the flux are correlated with the flux level itself in all energy band 
probed, indicating that lognormality is an instrinsic characteristic of PKS 2155-304. This results 
have already been shown for the flaring state of 2006 in the TeV band [4]. Others BE Lac type 
sources have present the same lognormal behavior, like BE Lacertae in the X-ray range [3] and 
Markarian 501 at TeV energies [16]. Interestingly non BE Lac objects such as the Seyfert 1 galaxy 
IRAS 13 244-3809 have revealed a similar behavior [17]. Seyfert 1 galaxies have different char¬ 
acteristics compare to BE Lac: radio quiet (no jet domination), with emission lines and more gas 
surrounding the central engine, indicating again that the lognormal behavior is disk related. 
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Figure 3: VED of PKS 2155-304. 


Moreover, the VED of PKS 2155-304 is showing an interesting feature with a similar variabil¬ 
ity increase in the X-ray and TeV ranges. Such a shape of the VED could help distinguish between 
different acceleration models in blazars as they won’t give the same variability signature. 
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